PROCEEDINGS 

0F SCIENCE 



Cold Nuclear Matter effects in T production in dAu 
collisions at RHIC 



(N 



A. Rakotozafindrabe 

IRFU/SPhN, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France 

E.G. Ferreiro 



DFP & IGFAE, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain 

3 F. Fleuret 

Laboratoire Leprince Ringuet, Ecole Polytechnique, CNRS/IN2P3, 91128 Palaiseau, France 

J.P. Lansberg 

^ IPNO, Universite Paris-Sud, CNRS/IN2P3, F-91406, Orsay, France 

Oh 

o | N. Matagne 

Universite de Mons, Service de Physique Nucleaire et Subnucleaire, Place du Pare 20, B-7000 
i Mons, Belgium 
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of the nuclear effects at play on top of the usual nuclear absorption, since the latter is expected to 
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the data leave some room for an additional fractional energy loss mechanism, recently revived in 
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1. Introduction 

Heavy-quarkonium (QQ) in-vacuum studies in proton-proton collisions allow to probe Quan- 
tum Chromodynamics (QCD) at short and long distances via their production mechanisms [1,2]. 
In-medium studies [2, 3] give access to Cold Nuclear Matter (CNM) effects in proton-nucleus col- 
lisions and to the understanding of QCD at high density and temperature in nucleus-nucleus colli- 
sions. STAR and PHENIX experiments [4] at RHIC provided the first measurements of the rapidity 
dependence of the inclusive T(bb) production in dAu collisions at ^/smn = 200 GeV. They are ex- 
pressed in terms of the nuclear modification factor Rj /Au of the T yield obtained in dAu collisions 
with respect to the superposition of the equivalent number of pp collisions at the same energy. De- 
spite the rather large uncertainties on /? J Au , we have shown [5] that it still gives qualitative insights 
into the nature of the CNM effects at play for T production in dAu collisions. 

At high energy (small Bjorken-jc), the nucleus in the rest frame of the incoming proton is 
Lorentz-contracted and the nucleon wave functions within this nucleus overlap. As a consequence, 
due to non-linear effects caused by the interactions of these overlapping nucleons, the nuclear Par- 
ton Distribution Functions (nPDF) deviate from those of free nucleons. Nucleons shadow [6, 7] 
each other and one expects the nPDFs to show lower values than for free nucleons. For ener- 
gies corresponding to 0.01 < xb < 0.3, some experimental data point [8] at an excess of partons 
compared to free nucleons, dubbed as anti-shadowing. 

In the following, we report on our recent study [5] where we have showed that usual nuclear 
modifications of gluon distribution in heavy ions - the shadowing and the anti-shadowing - as well 
as the possible break up (nuclear absorption) of the bb pair when it passes through the nucleus have 
a small effect: the resulting description of both the overall trend and magnitude of the T data in <iAu 
collisions appears rather unsatisfactory. Let us note that the mid rapidity measurement corresponds 
to the anti-shadowing region. 

Interestingly, we have pointed out that the modification of the T yield seen at backward ra- 
pidity corresponds to the intermediate Bjorken-x region, 0.35 < x# < 0.7, where a further nuclear 
suppression of the gluon distribution could be expected but was unobserved until now. Such an 
effect would be analogous to the one reported by the EMC collaboration [9] for the quarks. To 
date, there is no consensus [10] to explain this suppression, known as the (quark) EMC effect. It is 
still the object of intense investigations. This effect has not yet been clearly observed for gluons, 
even though it is allowed in some fits of gluon nPDFs. The gluon EMC suppression is usually 
overlooked whereas the shadowing of gluons is the subject of intense on-going discussions. Not 
much is known about gluons in this region and few data are sensitive to their distribution at xb 
larger than 0.3. The amount of the EMC suppression is actually basically unknown [1 1], except for 
a loose constrain set by the momentum sum rule. We explore the impact on Rj !Au that results from 
the use of the moderate and limiting cases currently offered by the nPDF parametrisations in this 
xb region. By comparing our results to the RHIC data, we can conclude that this provides a first 
and interesting hint for a strong gluon EMC effect, which might be stronger than the quark one. 

As regards the forward rapidity data, the discrepancy with the usual CNM effects cited above 
- shadowing and nuclear absorption - points at the presence of an additional CNM effect: a. frac- 
tional energy loss [12] proportional to the projectile parton energy and caused by medium-induced 
radiations associated to the quarkonium hadroproduction. This radiative energy loss arises when 
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the incoming partem and the outgoing coloured object radiate nearly coherently. This is possible 
when the colour-charge flow is subject to a scattering at a small angle in the nucleus rest frame. 
The bound on parton energy loss discussed previously in [13], which forbids energy loss to scale 
with energy, does not apply for such radiations. This fractional energy loss is thus probably at work 
at RHIC energies, contrary to the usual radiative energy loss [14]. Such an effect might induce an 
extra suppression of about 10 to 20% of the T yield in dAu at forward rapidity at RHIC. 

In practice, for the T production, we have used our well established Monte-Carlo framework 
JIN [15] - based on the probabilistic Glauber model and used to describe //\|/ production at RHIC 
- with the following ingredients: the partonic process for the bb production and the CNM effects. 
For now, the three T resonances are not resolved but are measured together. Since the nuclear 
absorption has to be small (as discussed in Section 4) and since the nPDF effects are very likely 
similar for these three states, we have safely considered them on the same footage. 

2. Partonic process for the bb production 

We have considered improved kinematics corresponding to a 2 — > 2 (g + g — > bb + g) partonic 
process for the T production. In earlier studies of CNM effects on T production [16], the bb pair has 
been assumed to be the result of a 2 — > 1 partonic process (i.e. g+g — >• bb). The presence of a final- 
state gluon introduces further degrees of freedom in the kinematics, allowing several (x\ , X2) for a 
given set (y, Pj) with the measured T rapidity and transverse momentum. Kinematics determines 
the physical phase space, but models are mandatory to compute the proper weighting of each 
kinematically allowed (x\,x%). This weight is simply the differential cross section at the partonic 
level times the gluon PDFs, i.e. g(xi,/UF )g(x2,HF)dG gg ^rc+gl dydPj dx\dx2. In the present status 
of our code, we are able to use the partonic differential cross section computed from any approach. 
For the present study, we have used the Colour-Singlet Model (CSM) at LO [17], which offers a 
good description of the direct T(15) (see Fig. la) and direct T(35) production at low Pj (where lies 
the bulk of the integrated cross-section). 

3. Gluon-monientuni distribution in nuclei 

To obtain the T yield in pA and AA collisions, a correction factor due to the nuclear modifica- 
tion of the gluon-momentum distribution has to be applied to the T yield obtained from the simple 
superposition of the equivalent number of pp collisions. This factor can be expressed in terms of 
the ratios Rf of the nPDF of a nucleon bound in a nucleus A to the free nucleon PDF. The numer- 
ical parametrisation of Rf (xb,Q 2 ) is given for each parton flavours. We have restricted our study 
to gluons since, at RHIC, T is essentially produced through the fusion of gluons [1,2]. As usually 
done, we label x\ (x%) the gluon momentum fraction in the proton/deuteron (nucleus). The nPDF 
spatial dependence has been included with a modification proportional to the local density [18]. 

To explore as widely as possible the impact of any nuclear modification of the gluon PDF, we 
have considered three different parametrisations: EKS98 [19], EPS08 [20] and nDSg [21] at LO. 
Yet, they span the current evaluation of the uncertainty on the gluon nPDF as provided by the newer 
sets, EPS09 [1 1], from a small to a very large (anti-)shadowing. The resulting /?J Au (see the green 
band in Fig. lb) has a much flatter rapidity dependence with respect to the data. The forward- 
y window, with (X2) ~ 0.008, corresponds to the shadowing region of the nPDF. However, the 
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Figure 1 



expected suppression is not strong enough to match the data. The mid-y window, with (jtj) ~ 0.05, 
corresponds to the start of the anti-shadowing region. More precise data is highly desirable to draw 
any conclsuion on any enhancement of the gluon nPDF in this region, the gluon anti-shadowing 
existence still being under debate. 

To what concerns the gluon EMC effect, we have used three of the EPS09 LO sets: one with 
a quark-like EMC gluon suppression, and the two limiting curves in the region 0.35 < xb < 0.7 
(see Fig. lc). They translate into the respective Rj lAu expectations at backward-j shown on Fig. Id. 
Despite the large experimental uncertainties, the comparison to the data already disfavours any 
nuclear PDF sets that enhance the gluon distribution in this xb region, while it favors a sizable 
gluon EMC suppression, analogous to the quark one or even stronger. 



4. Effective break-up cross-section 

The probability for the heavy-quark pair to survive the propagation through the nuclear medium 
is usually parametrised by an effective cross section a e ff. Due to its smaller size, the bb pair should 
suffer less break-up than the cc pair. Yet, the ratio of their size depends on the evolution stage of 
the heavy-quark pair. At the production time, this ratio is expected to be rti},/m c . When they are 
fully formed, it is rather ^4i"^ x as expected from their Bohr radii. The relevant timescale 



4 



CNM effects in T production in dAu collisions at RHIC 



A. Rakotozafindrabe 



for the pair evolution is the formation time. According to the uncertainty principle, it is related 
to the time needed - in their rest frame - to distinguish the energy levels of the IS and 2S states 
[22]: tf = ( M 2 2A/ ^2 \ ~ 0-4 fm for the T. For our purpose, tf has to be considered in the rest frame 
of the target (Au) nucleus. The formation time for different rapidities are given in Table 1: is 
significantly larger than the Au radius - except in the most backward region - implying that the bb 
pair is still in a pre-resonant state when traversing the nuclear matter. 

At forward and mid-y, this has two implications : cr J tf ~ {'^) 2 G J Jff ~ 0.1 following the 
early-time scaling m\,jm c \ and a e ff ought to be the same for the T(1S), T(2S) and T(3S) states, since 
they cannot be distinguished at the time they traverse the nucleus. On the contrary, at backward- 
y, we might expect different a e ff for these 3 states for y < — 1. However, the E772 experiment 
at Fermilab [23] has measured the T(1S) and T(2S+3S) separately at y/s = 38.8 GeV down to 
negative xp - with even smaller formation times - and it observed a similar suppression for the IS 
and the (2S + 3S) states. Such a result can only be understood if the absorption of the bb resonance 
is actually very small, preventing us to see any measurable difference between the 3 states. In 
the following, we will consider a range of G e ff from to 1 mb, where 1 mb is a conservative 
upper bound. In Fig. lb, this range translates into an uncertainty in Rj Au , which is obtained here 
using EKS98 as the nPDF parametrisation. Even in the strongest absoiption case, the resulting 
^-dependence of R^ Au is still not reconciled with the one of the data, which is much steeper at 
forward-^. 

Table 1 : T boost and formation time in the gold rest frame as a function of its rapidity at y/stw = 
200 GeV. 
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5. Conclusions and outlook 

We have reported on our recent investigations [5] on the T production in <iAu collisions at 
\/ s nn = 200 GeV, for which two RHIC experiments have provided the rapidity dependence. Due 
to the large scale set in by the T mass and the somewhat larger value of xi compared to the J j\f 
case, shadowing is found too weak to match the suppression seen at y > 0. In addition, the data at 
y ~ does not show any excess which would pin down anti-shadowing. Last, but not least, we have 
argued that the bottomonium survival probability to escape the Au nucleus should be quite large, 
compared to that of the charmonia at the same energy. Even in the strongest absorption case, the 
resulting y-dependence of Rj lAu cannot be reconciled with that of the data, which is more marked 
at forward-y. This discrepancy leaves some room for a fractional parton energy loss for forward 
angle in-medium T production, as recently revived in the literature [12] for the //\|/ production in 
pA collisions. 

In the most backward region, the suppression of the T yield may be the first hint of a gluon 
EMC suppression. Moreover, this gluon EMC effect might be stronger than the quark one which 
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still remains to be understood. Better precision measurements will be crucial to allow for a quanti- 
tative study of the gluon EMC effect, and will offer the opportunity to constrain the nuclear gluon 
distribution at large Bjorken-x. Along these lines, it has to be emphasised that a fixed-target exper- 
iment at the LHC, such as the project AFTER, would offer T yields in proton-nucleus collisions 
with nuclear targets three orders of magnitude larger than at RHIC [24, 25, 26]. Such precision 
studies of T production may in the future provide us with fundamental information on the internal 
dynamics of heavy nuclei such as those studied at RHIC and the LHC. 
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